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Abstract. We discuss the possible influence of fundamental QCD properties such 
as spontaneous chiral symmetry breaking and nucleon substructure on nuclear 
matter properties. We propose a chiral version of the relativistic a — to model in 
which the attractive background scalar field is associated with the chiral invariant 
field governing the radial fluctuations of the quark condensate. Nuclear matter 
stability is ensured once the scalar response of the nucleon depending on the 
quark confinement mechanism is properly incorporated. The needed parameters 
are estimated from lattice results and a satisfactory description of bulk properties 
follows, the only really free parameter being the uNN coupling constant. Pion 
loops can be also incorporated to obtain in a consistent way the finite density 
chiral susceptibilities. A good description of the asymmetry energy is obtained 
once the full rho meson exchange and Fock terms are included. 

1 Introduction 

A fundamental question of present day theoretical nuclear physics is the relationship between 
non perturbative QCD and the very rich structure of the nuclear many-body problem. A first 
trial beyond the standard non relativistic treatment based on effective forces (Skyrme or Gogny 
forces) is the relativistic mean field approach initiated by Walecka [T] . The next step should be 
to incorporate in this framework the effect associated with the most prominent properties of 
low energy QCD, i.e., the spontaneous breaking of chiral symmetry and the effect of hadron 
substructure and color confinement. In a bottom-to-top attitude this certainly constitutes a 
necessary basis to study matter under extreme density and temperature conditions. 



2 Chiral symmetry and confinement 

Basics on chiral symmetry. The SU{2)l SU{2)ii chiral symmetry of the QCD Lagrangian 
is certainly a crucial key for the understanding of many phenomena in low energy hadron 
physics. This symmetry originates from the fact that the QCD Lagrangian is almost invariant 
under the separate flavor SU (2) transformations of left-handed ql = (it/,, cIl) and right-handed 
Qb. = {ur, dft) light quark fields u and d. The explicit violation of chiral symmetry is governed 
by the quark mass vriq = {vriu -\- mdj/l < 10 MeV which is much smaller than typical hadron 
masses of order 1 GeV, indicating that the symmetry is excellent. It is however well established 
that the QCD vacuum does not possess the symmetry of the Lagrangian i.e., chiral symmetry 
is spontaneously broken (SCSB) as it is evidenced by a set of remarkable properties. The first 
one is the building-up of a chiral quark condensate : {qq) = {uu + dd)/2 which mixes, in the 
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broken vacuum, left-handed and right-handed quarks {{qq) — {qLqR + qRqL] /'^)- Another order 
parameter at the hadronic scale is the pion decay constant /tt — 94 MeV which is related to 
the quark condensate by the Gell-Mann-Oakes-Renner (GOR) relation : —2mg{qq)vac — "i^/^ 
valid to leading order in the current quark mass. It leads to a large negative value {qq)vac — 
-(240 MeV)'^ indicating a strong dynamical breaking of chiral symmetry. The second feature is 
the appearance of soft Goldstone bosons to be identified with the almost massless pions. Finally 
the chiral asymmetry of the vacuum associated with the condensation of quark-antiquark pairs 
is directly visible at the level of the the hadronic spectrum : there is no degeneracy between 
possible chiral partners such as p(770) - ai(1260) or 7r(140) - cr(400 - 1200). 

An illustration: the Nambu-Jona-Lasinio (NJL) model. Although chiral symmetry breaking and 
restoration can be studied on the lattice, the underlying physical mechanisms at the QCD level 
are not yet fully understood. It is thus useful to study this problem with models such as the very 
popular NJL model. In such a model the very complicated highly non perturbative multi-giuons 
exchanges between quarks are simulated by a very simple chiral invariant contact interaction 
with strength parameter Gi . In its simplest form the Lagrangian is : 

(-0^)" + (i^7^Ti/.)'] . (1) 

where ip = (ipuTipd) is an isodoublet of quark fields. The model also incorporates the fact that 
only low momentum quarks strongly interact in QCD. This is achieved by including a cutoff A 
of the order of 1 GeV in the momentum expansion of the quark fields. It is very easy to show 
that if the strength Gi is sufficiently large, the ground state wave function is of BCS type , 

= C exp ^— b^ps ^-p-s) I0o)j Sind made of interacting quark-antiquark pairs 

(created by the 6^ and operators acting on the bare perturbative vacuum \(f>o)), hence building 
up a quark condensate. As well-known such a BCS ground state is the vacuum of quasi-particles, 
the constituant quarks, with mass Mq solution of a gap equation : 

MQ=m-2G^{{qq))^m + 4N,G^ f 7^^, Ep^JIl^T?. (2) 

Jp<A \'^~) ^ 

It is thus clear that the existence of such a constituant mass of the order of 350 MeV is intimately 
related to the quark condensate {{qq)) = {{tpu'ipu + '4'dfpd))/'2 which plays the role of an order 
parameter associated with the broken vacuum. 

The mesons can be obtained using standard RPA exactly as in nuclear physics, as collective 
quark-antiquark excitations. This concerns in first place the pion (generated in the ij^j^rip 
channel) which appears as a Goldstone boson and its chiral partner the scalar-isoscalar sigma 
meson {ij;ip channel) which comes out with a mass TOq- ~ 2Mq. This model insisting on chiral 
symmetry, once taken in a more sophisticated form, yields a rather good phenomenology, espe- 
cially in the meson sector, but has an important weakness since it does not incorporate color 
confinement [2]. 

Interplay between chiral symmetry and confinement. When hadronic matter is heated and com- 
pressed, initially confined quarks and gluons start to percolate between the hadrons to be 
finally liberated. This picture is supported by lattice simulations showing that strongly in- 
teracting matter exhibits a sudden change in thermodynamic quantities (constituting a true 
phase transition or a rapid cross-over) within a narrow temperature window around Tc = 170 
MeV. This transition is accompanied by a sharp dropping of the chiral condensate indicating 
chiral symmetry restoration. However why these two transitions seem to occur simultaneously 
is still a major open question. This not fully understood duality between chiral symmetry and 
confinement is also present at the level of the nucleon structure and consequently in nuclear 
physics. Some of the pictures or models for nucleon structure put the emphasis on color confine- 
ment. This is particularly true for the MIT bag model where quarks with current quark mass, 
rriq ~ 5 — 10 MeV, move freely in a cavity (bubble of perturbative vaccuum) generated by the 
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Fig. 1. Left panel: Mexican hat effective potential. Right panel: Chiral circle corresponding to the bottom 
of the mexican hat potential 



confinement dynamics. A totally different view is provided by the constituant quark picture 
where the nucleon is made of three "big" quarks {i.e., of the NJL type) getting their mass 
from the chiral condensate of the non perturbative QCD vacuum. A plausible hybrid picture 
mixing both aspects seems to be supported by lattice calculations [S] . The nucleon looks like a 
Y shaped color string with three constituant quarks attached at the end-points of the Y. The 
introduction of a coupling between quarks and mesons allows to build the pion cloud which 
plays an extremely important role to account for the mass and the chiral properties of the 
nucleon (such as the sigma term or chiral susceptibilities defined below). Finally possible con- 
figurations where a string develops between one quark and one diquark in the color antitriplet 
channel, might also be present in the nucleon. 



3 Chiral effective theory: nucleon structure and nuclear physics 

How to build an effective chiral theory? NJL model as an example. For applications in hadronic 
or nuclear physics, one usually tries to formulate low energy QCD directly in terms of physical 
degrees of freedom, namely the hadrons. Mathematically one has to integrate out quarks and 
gluons in favor of hadrons. This is of course a formidable task that we will illustrate with the 
example of the NJL model. The appropriate technics is the path integral formalism where the 
partition function (vacuum persistence amplitude) is expressed as a functionnal integral over 
the quark fields: Z = J DipDip exp i^id'^xCjsijL)- The idea is to make a change of variable 
in the functionnal integral a = g iptp — 9 i'tpj^rij) {g is a coupling constant to be fixed) 
and to integrate out the quarks in the Dirac sea. One thus obtains a very complicated highly 
non local effective Lagrangian (the fermion determinant) depending on the mesonic fields, 
(T, TT, .... An approximative form can be obtained using a derivative or a loop expansion. To one 
loop, one gets: Leff — (1/2) {d^crd/j^a + d^n ■ dfj^ir) — [/(tr, tt). The effective potential U{(T,n) 
has the typical mexican hat shape associated with a broken (chiral) symmetry (see fig. 1). 
The degenerate vacuum corresponds to the chiral circle, i.e., the bottom of the mexican hat 
(fig.l). The pseudoscalar field (pion) represents the phase fluctuation of the condensate and it is 
massless since it does not cost energy to move on the bottom. The massive scalar mode (sigma) 
has a mass related to the curvature of the potential and represents the amplitude fluctuation 
of the condensate. The vacuum point with isospin zero is the point where tt = and where a 
has an expectation value identified with the pion decay constant /t^. 

Effective chiral theory and the background scalar field in nuclear matter. More generally we 
take the point of view that the effective theory has to be formulated first in term of the fields 
associated with the fiuctuations of the condensate parametrized in a matrix form W = a + ir-n. 
The sigma and the pion are promoted to the rank of effective degrees of freedom. An alternative 
and very convenient formulation of the resulting sigma model is obtained by going from cartesian 
to polar coordinates (right panel of fig.l), i.e., going from a linear to a non linear representation. 



4 



Will be inserted by the editor 



according to : 

W ^<T + IT -TT ^ SU = iU + s)exp(^^:^^ . (3) 

The new pion field corresponds to an orthoradial soft mode which is automatically massless 
(in the absence of explicit CSB) since it is associated with rotations on the chiral circle without 
cost of energy. The new sigma meson field 5', which is a chiral invariant, describes a radial mode 
associated with the fluctuations of the "chiral radius" It can be associated with the ordinary 
sigma meson which gets a very large width from its strong decay into two pions. Since it has 
derivative couplings, it decouples from low energy pions whose dynamics is described by chiral 
perturbation theory. With increasing density, its fluctuations s — S — fy^ are associated with the 
shrinking of the chiral circle and it governs the evolution of the nucleon mass. In a recent paper 
[4] we have proposed to identify this chiral invariant s field with the sigma meson of nuclear 
physics and relativistic theories of the Walecka type, or, said differently, with the background 
attractive scalar field at the origin of the nuclear binding. This also gives a plausible solution 
to the long-termed problem of the chiral status of Walecka theories. 

Tests of the effective theory with a chiral invariant scalar field. Once the appropriate couplings 
of the chiral fields to the baryons are introduced one can build an effective lagrangian to 
describe nuclear matter. Vector mesons {uj and p) must be also included to get the needed 
short range repulsion and asymmetry properties. At the Hartree level, the pion and the rho do 
not contribute for symmetric nuclear matter whose energy density written as a function of the 
order parameter s = (s) is : 

^p(^) — + M^(s) is the energy of an effective nucleon with the effective mass M^{s) = 
Mn + gs s. gs is the scalar coupling constant of the model; in the pure linear sigma model it is 
gs = Mn //,r ■ The effective potential U (cr, tt) when reexpressed in term of the new representation 
has the form : 

s is obtained by minimization of the energy density and is given at low density by : s ~ 
— {gs/'m'i) Ps- Its negative value is at the origin of the binding but the presence of the s'^ term 
(tadpole) has very important consequences. This tadpole is at the origin of the chiral dropping 
of the sigma mass Z\m* ~ —{2> gs /2f.,^)ps (a ~ 30% effect at po), and generates an attractive 
three-body force which makes nuclear matter collapse and destroys the Walecka saturation 
mechanism. Hence the chiral theory does not pass the nuclear matter stability test. 

This failure, which is in fact a long-standing problem [516] . is maybe not so surprising since 
the theory, as it is, also fails to describe some nucleon structure aspects as discussed below. 
The nucleon mass, as well as other intrinsic properties of the nucleon (sigma term, chiral 
susceptibilities), are QCD quantities which are in principle obtainable from lattice simulations. 
The problem is that lattice calculations of this kind are not feasible for quark masses smaller 
than 50 MeV, or equivalently pion mass smaller than 400 MeV, using the GOR relation. Hence 
one needs a technics to extrapolate the lattice data to the physical region. The difficulty of the 
extrapolation is linked to the non analytical behaviour of the nucleon mass as a function of mq 
(or equivalently m^) which comes from the pion cloud contribution. The idea of Thomas et al 
[7] was to separate the pion cloud self-energy, ST,{mT^^ A), from the rest of the nucleon mass 
and to calculate it in a chiral model with one adjustable cutoff parameter A. They expanded 
the remaining part in terms of as follows : 

MNiml) = flo + 02 774 + a4m'^ + Z'^(m^,yl) . (5) 

The best fit value of the parameter 04 shows little sensitivity to the shape of the form factor, 
with a value 04 2± —0.5 GeV~^ while 02 ~ 1.5 GeV~^ (see ref. [7]). The small value of 04 reflects 
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the fact that the non pionic contribution to the nucleon mass is almost linear in {i.e., in 
mq). Taking successive derivatives of Mjv with respect to {i.e., to niq), it is possible to 
obtain some fundamental chiral properties of the nucleon, namely the pion-nucleon sigma term 
and the scalar susceptibility of the nucleon. The non pionic pieces of these quantities are given 
by: 

cTjy — ~m^^ =a2TO^ + 2a4m^ ~29MeV. (6) 



o / non—pion /n \ 

dniq /4 dml \ ml j 



Xns H- -2 ,4 A^2 —2 =^^4a4. (7) 



In the above equations the first equalities correspond to the definitions, the second equalities 
make use of the GOR relation and the last ones come from the lattice QCD analysis. With 
typical cutoff used in this analysis, yl ~ 1 GeV, which yields crj^^ ~ 20 MeV, the total value of 
the sigma term is — 50 MeV, a quite satisfactory result in view of the most recent analysis. 
It is interesting to compare what comes out from the lattice approach with our chiral effective 
model. At this stage the only non pionic contribution to the nucleon mass comes from the scalar 
field, or more microscopically the nucleon mass entirely comes from the chiral condensate since 
the nucleon is just made of three constituant quarks with mass Mq = g{S)vac = s/tt — 350 
MeV. Hence the results for the non pionic sigma term and scalar susceptibility are identical to 
those of the linear sigma model : 

The identification of with cr^'' of our model fixes the sigma mass to a value iria- = 800 

MeV, close to the one ~ 750 MeV that we have used in a previous article ■ As it is the ratio 
gs/mi which is thus determined this value of ma- is associated with the coupling constant of 
the linear sigma model gs = / f-n = 10. Lowering gs reduces mo-. Similarly, the comparison 

of 

xivs '^'^^ ^^"^ lattice expression provides a model value for 04. The numerical result is 
a'^^ = — 3.4GeV~'^ while the value found in the expansion is only — 0.5GeV~'^. 



Nucleon structure effects and confinement mechanism. The net conclusion of the above discus- 
sion is that the model as such fails to pass the QCD test. In fact this is to be expected and even 
gratifying because it also fails the nuclear physics test. We will see that these two important 
failures may have a common origin. Indeed an important effect is missing, namely the scalar 
response of the nucleon, kns = B'^Mm /ds^ , to the scalar nuclear field, which is the basis of 
the quark- meson couphng model, (QMC), introduced in ref. The physical reason is very 
easy to understand: the nucleons are quite large composite systems of quarks and gluons and 
they should respond to the nuclear environment, i.e., to the background nuclear scalar fields. 
This response originates from the quark wave function modification in the nuclear field and 
will oviously depend on the confinement mechanism. This confinement effect is expected to 
generate a positive scalar response knSi it opposes an increase of the scalar field, a feature 
confirmed by the lattice analysis (see below) . This polarization of the nucleon is accounted for 
by the phenomenological introduction of the scalar nucleon response, K^rs, in the nucleon mass 
evolution as follows : 

Mn{s) = Mn + 55 s + - kats (9) 

This constitutes the only change in the expression of the energy density (eq. H]) but this has 
numerous consequences. The effective scalar coupling constant drops with increasing density 
but the sigma mass gets stabilized : 



(10) 
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Fig. 2. On the left panel the btndmg energy of nuclear matter with — 8.0, rric, = 850 MeV and 
C — 0.985 is shown. The full line corresponds to the full result, the dotted line represents the binding 
energy without the Fock and correlation energies and the dot-dahed line corresponds to the contribution 
of the Fock terms. The decreasing dotted line (always negative) represents the correlation energy. The 
right panel shows the density evolution of the QCD susceptibilities (normalized to the vacuum value 
of the pseudoscalar one). Dashed curve: pseudoscalar susceptibility. Full curve: scalar susceptibility. 
Dotted curve: nuclear contribution to the scalar susceptibility. Dot-dashed curve: pion loop contribution 
to the scalar susceptibility. 
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Fig. 3. Results of the RHF calculation with VDM parameters for the rho meson. Left panel: bin 
energy of symmetric nuclear matter. Right panel: density dependence of the asymmetry energy param- 
eter. 



The non-pionic contribution to the nuclcon susccptibiUty is modified, as well 



We see that the effect of confinement (kns) is to compensate the pure scalar term. Again 
comparing with the lattice expression one gets a model value for the 04 parameter : 

04 = - — (3 - 2C). (12) 
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where C is the dimensionless parameter C — (/^/2A/) kns- Numerically a4 = —0.5 GeV^^ 
gives C = +1.25, implying a large cancellation. As discussed below we will see that such a sig- 
nificant scalar response will generate other repulsive forces which restore the saturation mech- 
anism. At this point it is important to come again to the underlying physical picture implying 
that the nucleon mass originates both from the coupling to condensate and from confinement. 
In the original formulation of the quark coupling model, nuclear matter is represented as a 
collection of (MIT) bags seen as bubbles of perturbative vacuum in which quarks are confined. 
Thus in such a picture the mesons do not move and do not couple to quarks as in the true non 
perturbative QCD vacuum. Consequently the bag picture is at best an effective realisation of 
confinement which must not to be taken too literally. Indeed, QCD lattice simulations strongly 
suggest that a more realistic picture is closer to a y shaped color string (confinement aspect) 
attached to quarks. Outside this relatively thin string one has the ordinary non perturbative 
QCD vacuum possesing a chiral condensate from which the quarks get their constituant mass. 



4 Results for nuclear matter and discussion 

Stability and chiral properties of nuclear matter. In our first work devoted to the stability of 
nuclear matter |^ we limited ourselves to the Hartree approximation (eq. [3]) . The free sigma 
mass and the C parameter characterizing the nucleonic scalar response were taken as adjustable 
parameters. The fact that the best fit came with values very close to what we deduced afterwards 
from lattice is certainly satisfactory. One problem was the too large incompressibility K but 
is is nevertheless possible to improve our model description by adding an extra term in the 
expression of the effective nucleon mass in such a way that the nucleonic response vanishes at 
full chiral restoration (s = — /,r)- In this case the set of parameters : g^j = 6.8, rric — 750 MeV 
and C — I, leads to correct saturation properties, with an incompressibility value K — 270 
MeV. Another remarkable finding, according to the previous discussion, was the stability of the 
sigma mass which dropped by only 100 MeV at 2.5po- Ignoring the repusive tadpole associated 
with the nucleonic response would have given an almost vanishing sigma mass at this density. 

The next step jlO) has been to introduce pion loops which are necessary if we wish to address 
the question of the chiral properties of nuclear matter (namely the in-medium quark condensate 
and chiral susceptibilities) in a way which is consistent with nuclear matter stability since pion 
loops affect the energy. They do not enter at the Hartree mean field level but contribute through 
the Fock term and through the correlation energy. This energy can be calculated within a RPA 
ring approximation according to : 



i7L(w, q) is the full (RPA) longitudinal spin-isospin polarization propagator in the pionic chan- 
nel. This calculation actually includes iterated pion exchange and also the part of the NN 
potential from the two-pion exchange with A excitation. VL('i',q) is nothing but the (non 
static) pion exchange potential corrected by short-range components via the Landau-Migdal 
parameters, g^jvi s'na^ 9'aa- We have taken their values from a systematic survey of the data 
on spin-isospin physics by Ichimura et al. [11] : = 0.7, = 0.3, g'^^ = 0.5. Beside 

pion exchange we have also introduced the transverse channel (Vr), dominated by p exchange 
together with the short-range component. All the other needed ingredients are taken from the 
pion-nucleon and pion-nuclcus phcnomcnolgy. The scalar coupling constant is the one of the 
sigma model gs = Mn / f^^ — 10 and for the sigma mass we have followed the lattice indica- 
tions, allowing a small readjustment around the lattice value, which is m„ = 800 MeV. We 
have found a better fit with = 850 Me which corresponds to cr"°"-P»°" = 26 MeV. The 
omega mass is known and the loNN coupling constant is the only really free parameter. For the 
nucleon scalar response we have followed the indications of the lattice data but not strictly in 
view of the uncertainties attached to the higher derivatives. The value which fits the saturation 




([FL(c^,q)77L(t^,q;A)] + 2 [Vt{uj,(i) nT{io,q; 



A)]) 



(13) 
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properties is found to be C ~ 1, not far from the lattice value C = 1.25. With these inputs 
we have obtained a satisfactory description of the nuclear binding. The binding energy per 
particle is shown in fig. 2, with its different components. Taking the successive derivatives of 
the equation of state (the grand potential uj — e — ^p) with respect to mj, it is possible to 
obtain the pseudo-scalar susceptibility (directly related to the quark condensate) and the scalar 
susceptibilities of nuclear matter |12| . Since they are associated with chiral partners (pion and 
sigma), they should be identical at chiral restoration. The numerical calculation [lOj on the 
right panel of fig. 2 indeed displays a convergence effect. 

Relativistic Hartree-Fock (RHF). The description of asymmetric nuclear matter certainly ne- 
cessitates the full rho exchange already at the Hartree level through its time component. We 
thus include it using the following interacting lagrangian : 

Cp = ~gp paf. ^Tra^ - d,pa^ ^a^^'Ta^ + \ ml papp^, - i GTGaM. (14) 

One specific comment is in order for the tensor coupling of vector mesons. The pure vector 
dominance picture (VDM) implies the identification of Kp with the anomalous part of the 
isovector magnetic moment of the nucleon, i.e., Kp — 3.7. However pion-nucleon scattering 
data suggest Kp = 6.6 (strong rho scenario). At the Hartree level we still have our previous 
picture. Symmetric and asymmetric nuclear matter are seen as an assembly of nucleons, i.e., 
of Y shaped color strings with massive constituant quarks at the end getting their mass from 
the chiral condensate. This nucleons move in self-consistent scalar (s) and vector background 
fields (w^jPg). However to get a quantitative description, all the Fock terms, including not 
only pion and rho exchange but also sigma and omega exchange, are needed. These exchange 
terms are mediated by the propagation of the fluctuations of the previous meson fields. In 
particular it is possible to show that the scalar fluctuation, i.e., the fluctuation of the chiral 
condensate, propagates with an in-medium modified sigma mass m*^ = ^"(s) + kj^s PS- The 
results of preliminary RHF calculation [13] with rho coupling given by pure VDM {gp = 2.65, 
Kp = 3.7) are shown on fig. 3 and reproduce the saturation properties, g^j = 6.35 is not very 
far from the quark model value g^j = Sgp and C ~ 1.75 is larger than the lattice value. The 
resulting incompressibility K — 242 MeV and asymmetry energy parameter as = 30.2 MeV 
are obviously satisfactory. It is also possible to reproduce saturation properties in the strong 
rho scenario {Kp = 6.6) but the asymmetry energy becomes too large : as = 45 MeV. These 
results are certainly very encouraging but the ultimate test of this chiral theory will be the 
reproduction of finite nuclei (in particular the spin-orbit splitting) and the consequences for 
neutron star properties. 
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